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I INTRODUCTION 

The inverse Compton <IC) process occurs when low frequency 
photons are scattered to higher frequencies by relativistic 
electrons of Lorenz factor y , with an energy boost of ~y . The lobes 
and jets of radio galaxies contain huge reservoirs of synchrotron 
emitting electron with ~ 1 0 2 - 10 4 , It is thus expected that they 
will scatter microwave background photons with v ~ 1 0 1 * Hz into the 
X-ray band where they would be seen by satellite bourne X-ray 
detectors. Despite the obligatory nature of this process, there has 
not been a single confirmed detection of X-ray IC emission from 
extended extragalactic radio structures. The most likely case prior 
to the launch of the Einstein Observatory, a report of X-ray emission 
from the outer lobes of Centaurus A from the Ariel V satellite (Cooke 
et al. 1978), was shown by SAS-3 researchers to be due to confusion 
by an unrelated X-ray source (Marshall and Clark 1981). Although a 
few jets and lobes appear to have X-ray counterparts in Einstein 
images (e. g. , Feigelson et al. 1981; Harris et al. 1984; Burns et al. 
1985; see reviews by Feigelson 1983 and 1985), they can not be 
confidently attributed to IC processes. The bright X-ray jet knots 
of Cen A and Vir A, in particular, are most likely due to direct 
synchrotron from extremely energetic electrons (Schreier et al, 

1 982; Burns et al. 1 983). 

This general failure to detect IC X-rays appears to be due to a 
conspiracy of numerical factors. Given the estimated magnetic field 
strengths of B eq ~ 1 0" 6 - 10"^ G and the microwave background photon 
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density U ph ~ 10 erg cm , typical radio components with S v < 100 
Jy at 100 MHz will produce X-ray flux densities < 0.001 - 0.0001 yJy, 
which is below the sensitivity limits of current X-ray 
instrumentation (Harris and Grindlay 1979). Only unusually strong 
radio sources with unusually weak magnetic fields should be 
detectable. Current upper limits on IC X-rays generally provide only 
uninteresting limits on the field strength (Feigelson and Berg 1 9 6 3 > • 
with a few exceptions such as the Cen A outer lobes limit (Marshall 
and Clark 1981). A single confirmed detection of IC would be highly 
significant, as it would determine whether the many assumptions 
implied by the use of equipartition between the magnetic field and 
electron pressures (e. g. Miley 1 980) are in fact valid. These 
include assumptions regarding the pressure of relativistic protons, 
the filling factor, the orientation of the magnetic field, and the 
source geometry. 

A number of researchers (Turland 1 9 7 5 i Andernach et al. 1 979; 

Dennison 1 980) have noted that the radio halo of Virgo A = 3C274 = 

M87 provides an excellent opportunity to detect IC X-rays because of 

its high flux (S^ = 1360 Jy at 101 MHz) and low equipartition field 

strength (B = 4-8 PG). The predicted IC flux would be of order 
eq 

10 of the thermal emission from the dense 'cooling flow' of 
intracluster gas that is responsible for most of the X-ray emission 
of M87 (e. g. Fabricant et al. 1 980; Stewart et al. 1984). It might 
be resolved from the thermal component because of the distinctive 
morphology of the radio halo: two misaligned lobes 2' to 6* east and 

south-west of the nucleus (Andernach et al. 1979; Kotanyi 1980). 
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When the Einstein High Resolution Imager (HRI) obtained a deep 
exposure of M87, an X-ray component with just this morphology was 
seen superposed on the large-scale thermal component (Schreier et al. 
*' 962 , henceforth SGF). The observed excess of ~0. 02 HRI cts s ^ was 
within a factor of two of the predicted IC level, and thus represents 
the best case to date for X-ray IC emission from an extended radio 
structure, 

The purpose of this study is to investigaate in more detail the 
associated X-ray and radio emission in the Virgo A halo discovered by 
SGF. Improved Einstein HRI data and new radio maps obtained with the 
Very Large Array are described (section II), and the relation between 
the X-ray and radio structures is carefully examined (section III). 
Section IV presents several possible explanations for the X-ray 
emission. The IC model is found to be viable only if the magnetic 
field is variable and substantially weaker than the equipartition 
value. The principal alternative is excess thermal X-rays due to 
compression of the intracluster medium by the radio lobe. In either 
case, the association of such prominent radio and X-ray structures is 
unique among known radio galaxies. 

II Observations and Data Analysis 
(a) X-ray Observations 

The study by SGF was based on a 75, 600 s exposure of M87 
obtained on July 5 1979 with the Einstein HRI. The instrument is 
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described by Giacconi et al. (1 979). We have added to this image a 
second exposure of duration 51 , 200 s obtained on December 9 1980. 

Both data sets were reprocessed with improved aspect solutions (Rev. 

1 . 0 > i and data acquired when the star trackers were out of LOCK mode 
were rejected. The resultant merged image contains 7*4* H5 s 
effective exposure and should be free of aspect errors greater than 
1-2 arc seconds. The HRI is not subject to any other temporal or 
spatial anomaly; detector background is ignored because it is nearly 
uniform and will be removed in the analysis of asymmetric emission. 

As described by SGF, the image shows four X-ray components; a 
cusp or point source associated with the active nucleus; a bright 
region from the radio/optical synchrotron jet 10”-20” west of the 
nucleus; diffuse emission extending beyond the 35* HRI field 
corresponding to the dense thermal gas; and the asymmetric excess 
around 3' - 5' east and south-west of the nucleus associated with the 

radio halo. Figure la shows a contour map of the merged HRI image 
after smoothing with a 20” FWHH Gaussian. The slightly trangular 
contours are the only sign of the asymmetric component in this 
representation. 

The main task in further revealing the structure of the 
asymmetric component associated with the radio halo is to remove the 
thermal gas component. Though in principle one could subtract a 
theoretical model of the gaseous emission based on other X-ray data 
and a model of the gravitational potential, in practice researchers 
have assumed the X-rays are entirely thermal in origin and have used 
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them to derive the mass distribution <e. g. Fabricant et al. 1980; 
Stewart et al. 1 98 H > . Since there is no adequate independent 
knowledge of the mass distribution, we have made the a priori 
assumption that the thermal component contributes all of the 
spherically symmetric HR1 counts, with the additional condition that 
its radial distribution monotonically decrease with distance from the 
nucleus. Any residual asymmetric emission is then compared to the 
radio structure. There is no guarantee that the gaseous component is 
truly symmetric (indeed, we attribute at least some to the asymmetry 
to it; see Appendix A), or that the symmetric component does not 
include some nori-thermal emission* Despite those limitations, the 
separation between symmetric and asymmetric components suffers little 
operational ambiguity. 

To obtain a map of the asymmetric X-ray component, the HRI image 
was binned into 4" pixels and divided into azimuthal octants and 
differentially spaced annuli centered on the nucleus. The inner ^0" 
with the jet and steeply rising core region were ignored. The octant 
in each annulus with the lowest surface brightness was chosen to 
represent the symmetric component at that radius. The lowest octants 
generally occurred with P. A. 0 - 90d. The resulting surface 
brightness distribution follows closely the lowest curves shown in 
Figure 3 of SGF. Using parabolic interpolation between the annuli, 
the symmetric surface brightness was evaluated and subtracted from 
each pixel in the map. Figure 1b shows the residual map after 
smoothing with a 15" FWHM Gaussian. 
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conversion from HRI cts s -1 to physical units is quite uncertain 

given the lack of energy discrimination in the HRI and the unknown 

spectral characteristics of the asymmetric component. We have 

adopted a spectrum expected if the emission is truly IC: a power law 

with energy index a = 1 . 0 (taken from the radio spectrum around 

100-400 MH 2 ; Andernach et al. 1979) and a galactic absorption of N„ = 

H 

4 x 1 0 20 cm” 2 . The conversion factors between 1 HRI count/1000 
seconds and physical units are then (Harris 1984) 1 ct/1000 s = 0.016 
yjy (evaluated at 1 keV) and 1 ct/1 000 s = 2. 9 x 1 0 39 erg s” 1 
(assuming a distance of 1 6 Mpc). If the emission process is thermal 
bremsstrahlung from gas with cosmic abundances at T = 2 x 10^ K, 
which is the estimated temperature for gas several arc minutes from 
the nucleus (Stewart et al. 1984), the flux density and luminosity 
conversion factors would be 5.7 and 1.2 times higher respectively. 

The flux density is thus extremely sensitive to the unknown spectral 
index of the asymmetrical component, though the broad band flux and 
luminosity is not. Ue thus caution that, while the residual map 
shown in Figure 1b is an accurate representation of the observed 
asymmetric X-ray emission, the associated flux density map might 
differ significantly due to spectral variations. 

(b) Radio Observations 

Given the excellent X-ray residual map with resolution of 15" or 
better, we sought to obtain a radio map of the Vir A halo with 
similar resolution. This is extremely difficult because the halo has 
a very low surface brightness and surrounds a very bright core. (The 
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cope component is itself complex wi*,h a VLBI source at the nucleus, 
the one-sided jet, end two inner lobes. ) There is an unavoidable 
tradeoff between resolution and dynamic range; dynamic ranges of 10 3 
- 10^ are needed at sub-arcminute resolution. The effort is further 
complicated by the large angular size of the halo (around 15') which 
can be overresolved by widely spaced interferometers and can exceed 
the primary beam widths of some telescopes. Observations of the halo 
are best done at low frequencies where its flux density exceeds that 
of the core component. 

There have been three efforts to map the two-dimensional 
structure of the halo, which was discovered by Mills (1952) in an 
early interferometric experiment at 101 MHz. A crude interferometric 
map with 3' x4’ resolution was made by Cameron (1971) at 408 MHz with 
the Molonglo Cross. Andernach et al. (1 979) obtained a 4. 8 GHz map 
with 2.6' resolution and a 1 0. 7 GHz map with 75'' resolution with the 
single Bonn 100-m telescope. The latter shows few details because of 
the low halo flux at this high frequency. Lastly, Kotanyi (I960) 
made a 610 MHz map with 4. 3' x0. 9* resolution with the Westerbork 
interferometer. The elongated beam is due to the absence of 
north-south baselines. Polarization maps were also obtained from the 
Bonn and WRST data. 

The NRA0 Very Large Array (VLA) in its low resolution 
configurations should provide an improved view of the halo. The 
instrument is described by Thompson et al. (1980). Three data sets 
were available for investigation: (1)A full synthesis at 1 . 4 GHz, 
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interspersed with snapshots at 5 GHz, obtained on January 30 1982 
with a hybrid C/A configuration (north arm extended); (2) full 
syntheses at 1,4 and 5 GHz obtained on April 6 1 9 8 *4 with the 
standard C configuration; and (3) a series of snapshots spanning 6 
hours with 1 hour net exposure at 1 . >l GHz and 5 GHz respectively 
obtained on June 2 1983 in the D configuration, The second data set 
was a byproduct of an intercontinental VLBI observation of the Vir A 
nucleus (see Reid et al, 1982 for an earlier epoch of this 
experiment). The third data set was part of an effort to map the 
calibrator sources of Baars et al. 0 977). and was kindly made 
available by Dr. R. Perley. The antenna phases and gains were 
calibrated for each data set using standard VLA procedures. Flux 
densities were bootstrapped to 3C286, Maps were made for each 
database using a variety of tapering and weighting functions, after 
which there were CLEANed, The visibility data were then 
self-calibrated to a model based on the clean components, 

The difficulties of mapping the halo at high resolution were 
rapidly encountered. The 1.H GHz C/A hybrid and C configuration 
maps, with synthesized beamwidths between 5" and 30" depending on the 
tapering function, could not be adequately cleaned or 
self-calibrated. Dynamic ranges reached 1000:1 (Peak flux 
density: RMS noise) but the low halo surface brightness at these high 
resolutions requires dynamic ranges several times higher to see 
details. Furthermore, only ~70% of the 200 Jy in the source appears 
in the cleaned map because much of the halo flux occurs on very low 
spatial frequencies. It is resolved out in the C configuration with 


ft minimum baseline of 31 5 wavelengths, Insertion of a large 
zero-spacing flux does not improve map quality, probably because of 
the non-uniform structure of the large scale emission. 

The only acceptable 1 . 4 GHz maps were obtained with the D 
configuration with a minimum baseline of 200 wavelengths. The data 
were subjected to four pases of cleaning and self-calibration, in 
which phases only were adjusted. The recovered flux density is 197 
Jy which is consistent with the total 209 + /- 9 Jy in the source at 
1.465 GHz (Bears et al. 1 977), Figures 2a and 2b show two 
representations of the resultant total intensity map, The first is 
subjected to a Gaussian taper at 1500 wavelengths and has a 
synthesized beam width of 72” FWHM, The second is untapered with 
'super-uniform' weighting in the UV plane (see Sramek 1 983) to give a 
high resolution of 39” FWHM. The dynamic ranges of the maps are 
1 300:1 and 3 M 00: 1 respectively. The polarization properties are 
shown in Figures 2a and 2b as electric field vectors proportional in 
length to the polarization fraction at each location. 

In addition to the 1.4 GHz maps, the 5 GHz data were also 
examined. The three-fold higher resolution, lower halo flux and 
smaller VLA primary beam (9* FWHM at 5 GHz) all conspire against 
retrieval of the more extended parts of the halo. Most the 5 GHz 
data sets failed to show more than 10% of the estimated 16 Jy in the 
halo (Andernach et al. 1979) due to the missing short baselines, 
nonetheless, we show the best obtained map,, from the C configuration 
data tapered at 3000 wavelengths, in Figure 3. It has a 40” 
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synthesized beamwidth and shows a dynamic range of "SODOM, A 
significant primary beam attenuation correction factor has been 
applied. No polarization data is available beoause the array was 
phased for the VL.BI experiment. We emphasize that, beoause of the 
missing flux, the brightness of the features in Figure 3 should not 
be used for quantitative measures of flux or spectral Indices. 

Ill Results 

(a) X-ray and Radio Morphologies 

The asymmetrical X-ray emission appears on two spatial scales in 
Figure 1b; the bright emission within S' of the nucleus associated 
with the radio halo; and a low surface brightness excess seen out to 
"10* from the nucleus to the NW and SE. The latter Is seen in the 
contour map as a sprinkling of low level contours, and is quite 
apparent on displays at lower resolution. Unrelated to the radio 
source, this X-ray structure is probably associated with an 
eccentricity in optical isophotes seen in M87. Since this diverges 
from the principal topic of this paper, we discuss it briefly in 
Appendix A. 

The X-ray component associated with the radio halo might be 
described as a triangular region with the nucleus at one vertex and 
two adjacent arms oriented along P. A. 90° <4' long) and P. A. 225 ° < 5 ' 
long). The region between the arms is smoothly filled, with a 
surface brightness 50-80% that of the arms. The eastern aim is 



12 


* 


* 


* 


* 


brightest 2.0* from the nucleus, Its width st this point is 1.3' 

FVHM end 2. 5' FUZI. The arm dims end tapers to the east) ending at a 
resolved spot 1.0' FVHM centered H.1’ from the nucleus along P, A, 93*. 
The southwestern arm also peaks 1.0* from the nucleus, but has a 
somewhat different transverse brightness distribution, The north 
side of the structure is quite straight and sharp with 0. 3' HUHN, 
while the south side declines more gradually with 1,5' HVHM, The 
sharp edge of this arm is visually quite evident even in the original 
HRI image dominated by the symmetrical thermal component. 

A matter not fully resolved is whether the X-ray asymmetries 
extend directly into the nucleus. Pixel digitization and the steep 
gradient of the subtracted symmetric component hinder detailed 
analysis. Nonetheless, it appears that the eastern arm may extend 
directly into the central 20*' along P. A, 110°, while the southwestern 
arm appears to start about >10" from the nucleus. 

The radio structure shown in Figures 2a and 2b reveal 
considerably more detail than earlier maps of the halo, The total 
extent of the source is 13' elongated along approximately P. A. *»5°. 

But, unlike most double sources the interior features are generally 
not aligned along this axis. The principal features might be 
described as follows (in order of increasing distance from the 
nucleus): (a) an unresolved knot 2.1' west of the nucleus along P. A. 

265°; (b) an amorphous component to the east centered 3.0' from the 

nucleus along P. A, 85°, 2.5' FWHM in size; (c) a jet-like component 
extending 2* to 5* southwest of the nucleus along P. A. 15 s , with 


width 1.6' FWHM transverse of the axis; and (c) the large-scale 
envelope mentioned above. It is not clear these features are truly 
independent of each other. The higher resolution maps »t M and 5 
GHz (Figures 2b and 3) distinctly suggest that emitting material 
flows through the western f.tnot into the southwest feature as a 
continuous! knotty jet. The southwest jet fades smoothly into the 
large-scale envelope, though the latter is present on the northeast 
side of the source where no jet is seen. 

The polarization structure of the halo is even more unusual. 
Though ~H0% .of the unpolarized halo flux lies east and north of the 
nucleus, less than 5% of the polarized flu:< lies in these directions, 
Virtually all of the polarized emission is concentrated at the far 
end of the southwest jet. Its principal component is a ridge 4' in 
length oriented along P. A. 135 s , or transverse to the axis of the 
unpolarized emission. Yet, the direction of the magnetic field 
appears generally parallel to the unpolarized axis, which is contrary 
to the orientation generally seen in low luminosity radio jets 
(Bridle 1984). Some weak polarized emission is 3een near (but not 
coincident with) the west knot, again with a field orientation 
suggesting a flow from the knot into the southwest jet. All 
Interpretation of the magnetic field orientations is somewhat 
speculative, as the importance of Faraday rotation at 1 . 4 GHz has not 
been established. 

Xn Appendix B we discuss the halo radio morphology in light of 
current models of radio galaxy structure. We find the halo can 
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probably be classified as a wide-angle tail (WAT) source, Its 
structure might be attributed to a rotating nucleus, but more likely 
to bent Jets seen in projection. Further study, particularly 
multi-frequency polarization maps to eludicate the Faraday effects, 
is warranted. 

The descriptions given above indicate that the morphologies of 
the halo at the two wavebands are similar but not identical. The 
resemblances lie in the principal axes of emission, east and 
southwest of the nucleus, and the overall extent of the emission. 

Diff erences . include the brightness profiles along the principal axes, 
(e« g. , the X-ray emission peaks closer to the nucleus than the radio 
emission), and radio structures that do not appear in the X-ray map 
<e. g. the western knot and large-scale envelope to the north). The 
radio map resolution is not high enough to say whether or not radio 
emission is absent to the southeast of the nucleus where the 
prominent X-ray bridge appears. The only case of directly coincident 
features at the two wavebands are the east radio peak and the faint 
X-ray spot at the end of the eastern X-ray arm. 

(b) Quantitative Comparison of the X-ray and Radio Naps 

We first consider the halo structures as a whole. Within the 
lowest contour of Figures 1b and 1c, the HRI residual image shows 
about 10,300 counts or 138 cts/1000 s in a 19 square arc minute 
region. This value is accurate to about 10%, and excludes the inner 
HO" core region. This count rate of the asymmetrical component 
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comprises about 25% of the count rate seen in the original HRI image 
(Figure la). It is more than 3 times higher than the 40 cts/1000 s 
estimated by SGF from a less detailed analysis because they 
considered only a limited range of radii (2' to 5* ) and position 
angles (90° - 120° and 210° - 240°). Using the conversions given in 
section Ila, this count rate corresponds to a 0. 2-4 keV luminosity of 
4. 0 x 1 ergs s“* (accurate to 20% except for possible distance 
uncertainties) and a 1 keV flux density of 2. 2 pJy (which may be 
inaccurate by a factor of 5 due to spectral uncertainties). 

Within, the lowest contour of Figures 2a and 2b, the 1.4 GHz 

radio flux density is 197 Jy of which 58 Jy lies outride the central 

core region. The separation of the core from the halo is not 

clearcut; Andernach et al. (1 979) show than a decomposition of 120 

Jy in the core and 80 Jy in the halo give more consistent 0.1-10 GHz 

spectra for both components. Assuming the latter value, then the 

X-ray/radio flux density ratio log [S (1 keV) /S (1 . 4 GHz)] = -7.56 

corresponding to a spectral index 

a = -0. 92. 
rx 

As discussed in the previous section, the relative radio/X-ray 
emission is highly variable across the source. Figure 4 shows a 
superposition of the 39" resolution radio and X-ray residual maps 
with eight locations marked that show a range of different 
properties. Table 1 gives the measured radio flux density, 
polarization fraction, HRI counts/39" beam, estimated X-ray flux 
density, X-ray/radio flux density ratio and spectral index for each 


location. The results may be summarized as follows. The radio 
surface brightness ranges over a factor of 5 and the polarization 
fraction varies between 28% and <1% (three-sigma limit). The X-ray 
surface brightness ranges over a factor of >7. The X-ray/radio ratio 
varies by a factor >10 across the halo, with a spectral index of 
approximately a rx = -0.9 to -1.0. 

IV DISCUSSION 

Before considering the possible explanations for the associated, 
but not entirely coincident, X-ray and radio structures in Vir A, we 
briefly compare it to other sources. In both wavebands, the 
structure is remarkably bright. In X-rays, its HRI count rate is 5 
times greater than the jet in Vir A (SGF) , 20 times greater than the 
jet in Cen A (Feigelson et al. 1981), 70 times greater than the jet 
in 3C273 (Willingale 1981), and an order of magnitude brighter than 
possible X-ray counterparts of radio structures in other galaxies 
(Harris et al. 1 9 8 »i i Burns et al. 1 985). Its soft X-ray flux is 
sufficient to have been detected in the UHURU, Ariel V and HEA0-1 
all-sky surveys even without the presence of the dominant symmetric 
component. In radio, the Vir A halo is the third brightest 
extragalactic radio source at 178 MHz. Yet, if M87 lay at a distance 
of 500 Mpc (z = 0.10 - 0.15) rather than 15 Mpc, only 5-10 counts 
would have appeared in a long HRI exposure and the X-ray source would 
be only marginally detectable. At that distance, its 1.4 GHz radio 
flux density would be 80 mJy. Since there are large numbers of such 
weak distant radio galaxies, it seems quite likely that many similar 
X-ray structures are present at levels below that detectable with the 
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Einstein Observatory. 

(a) Inverse Compton (10) Emission 

As discussed in section I, the radio halo’ s relativistic 
electrons must scatter photons of the microwave background into the 
X-ray band. The ratio of IC to synchrotron flux density from this 
process, assuming a background blackbody spectrum with T = 2. 7° K and 
ignoring negligible redshift factors, is (Harris and Grindlay 1979) 

log(S /S ) = -16.24 + 4.70 a- a log(v /v ) - (ct + 1) log B (1) 

XV X IT 
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where v x is the observed X-ray frequency of 2.4 x 10 Hz (1 keV), 
is the observed radio frequency of 1.4 x 10 9 Hz, a is the spectral 

index of the (assumed) power law synchrotron spectrum (S ~ v a ) , and 

v 

B is the magnetic field strength in gauss. Application of this 

formula is complicated by the curvature of the Vir A radio spectrum, 

with spectral index ranging from a = 0.9 at 100 MHz to a = 2.3 at 10 

GHz (Andernach et al. 1979). We must estimate the effective 

spectral index between the observed frequency 1.4 GHz and the 

frequency emitted by electrons responsible for the IC scattering of 

microwave photons to 1 keV. These electrons have Lorenz factors y ~ 

(2 x 1 0 1 ? Hz/2 x 1 0 11 Hz) ~ 1000 and emit synchrotron around v ~ 80 
2 

(B/5 yG) (y/1000) MHz. Extrapolating the halo spectrum of Andernach 
et al. below 101 MHz, we obtain an effective spectral index between 
1.4 GHz and 80 MHz of a = 1.05. Using this value in equation (1) 
gives a predicted IC flux 
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log<S x /S v > = -9.08 - 2.05 log (B/5 y G) (2) 

The observed flux density ratio averaged over the entire halo is 

log (S^/S^) = -7.96 or 13 times this predicted level. It ranges 

between >6 and 70 times the predicted level at different location 

across the source (Table 1). This discrepancy is disturbing, but not 

conclusive, evidence against the IC model, The observed ratio S /S 

x v 

is unlikely to be low by more than a factor of 1.5 - 2. The 
conversion between HRI counts and flux density can only be ~50% 
higher than the value adopted in section Ila for any plausible power 
law plus absorption spectrum (Harris 1984). (The factor of 5 
uncertainty in conversion factor mentioned in Ila applies only if a 
thermal spectrum is permitted. ) The uncertainties in HRI counts 
rates or radio flux densities are probably below 20%. The 
discrepancy must therefore be associated with the predicted ratio in 
equations (1) and (2). The effective spectral index is accurate to 
+ 0.1, and the constants in equation (1) are accurate to 20% (Harris 
and Grindlay 1979). 

The only was to save the simple IC model is to invoke a magnetic 

field considerably below the average equipartition level of 5 yG. An 

average level of 1. 5 yG across the halo and, at the east X-ray peak 

where log (S /S ) is highest, a level as low as 0. 6 yG is needed to 
xv 

get the observed high X-ray flux. The equipartition level of 4 - 6 
yG was estimated by Andernach et al. (1 979) and Kotanyi (1980) from 
the overall properties of the halo. Although we do not have 
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localized spectral index maps or knowledge of the line-of -sight 
depths of the source at different locations, we can estimate the 
equipartition pressure in the west radio knot and southwest radio 
peak to be 13 yG. (A flux density of 1 Jy at 1J GHz in a 39"x39" 
beam is observed; a 39" line-of-sight depth and as 1.0 is assumed.) 
The magnetic field must thus be ~30 times below the equipartition 
level to account for the high X-ray flux at these locations. 

(b) Thermal Emission 

The principal alternative to the IC model is that the asymmetric 

X~ray component is due to thermal bremsstrahlung just like the 

symmetric component. Since thermal X-ray emissivity is strongly 

dependent on density and only weakly dependent on temperature, this 

requires a concentration of dense plasma at the locations of the 

X-ray structures east and southwest of the nucleus. Any local 

heating of ambient medium without simultaneous compression, due for 

instance to the leakage of relativistic electrons from the radio jets 

<e. g. Sofia 1 973; Scott et al. 1 980), is thus unlikely to account for 

41 -1 

the X-ray excess. The observed X-ray luminosity of 4 x 10 ergs s 

3 3 

is generated in a volume is ~5 x 10 kpc , assuming the line of sight 

depth of the region is similar to its north-south extent. The 

-2 -3 

required density is then ”1 x 10 cm for any reasonable 

6 7 0 

temperature between 5x10 and 5x10 K. The external medium in 

the vicinity of the asymmetric component has density ranging from 0. 3 

-4 -3 7 

-3x10 cm and temperature around 2 x 10 °K (Stewart et al. 

1984). Thus, if the asymmetrical emission is thermal its pressure 
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mu3t exceed that or the ambient medium by one or two orders of 
magnitude. Assuming an equipartition magnetic field, it will exceed 
the internal pressure of the radio emitting region by a similar 
factor. 

Given this substantial overpressure, the thermal model implies 
that the emission is transient. The gas will simultaneously cool due 
to line and bremsstrahlung emission, and expand into the surrounding 

q 

gas on timescales of ~1 x 10 years. The exact balance between these 
processes depends on its unknown temperature; cooling will proceed 
faster than . dissipation if the temperature is less than ~1 x 10 7 °K. 
Though it is not self-evident whether this gas is entrained within 
the radio emitting region or lies outside of it, the latter 
possibility seems more likely. If entrained, the gas must be 
inhomogeneously distributed within the radio Jets and would disrupt 
the field geometry due to its higher pressure. It would naturally 
occur outside the radio region if the jet* s outward substantially 
exceeds the ambient sound speed of 500 km s”* . The external medium 
would then be shocked and swept up by the jet much like a supernova 
remnant sweeps up the interstellar medium during its ' snowplow' 
phase. A major advantage of such a model is that precise spatial 
coincidence between the X-ray and radio emission is not expected. 

The fact that the X-ray component is brightest closer to the nucleus 
than the leading edges of the radio components may be due to the 
higher ambient gas density closer to the nucleus. 


(c) Other Models 
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The third possible emission mechanism for the asymmetric X-ray 
component, synchrotron emission, can be readily discounted. The Jet 
knots of Vir A and Cen A, where the X-rays are thought to be 
synchrotron, have flat radio spectra (a = o. 6) extending through the 
gigahertz band and (at least for Vir A) into the optical band. The 
Vir A halo has a rapidly steepening radio spectrum and no associated 
optical emission (Carter and Dixson 1978). Extrapolation of the 
radio spectrum into the X-ray band predicts flux densities orders of 
mangitudes below the observed levels. 

Another possibility is that the X-ray emission has no relation 

to Vir A = M87i but is a chance superposition of a foreground or 

background object. Huchra and Brodie (1984) have found several 

galaxies within 3* of the center of M87 with redshifts of 25800 km/s. 

They appear to be members of a group or poor cluster 25 times more 

distant than the Virgo cluster. Huchra and Brodie suggest that the 

intracluster medium of this cluster may be responsible for the 

asymmetric X-ray emission, in which case its soft X-ray luminosity is 
44 -1 

~2 x 10 erg s . We believe this explanation is improbable. 

First, this luminosity is very high, though not unprecedented, for 
poor clusters (Kriss et al. 1983). Such high luminosities only occur 
in cD-dominated poor clusters which have spherically symmetric X-ray 
distributions. Poor clusters with lumpy and asymmetric gas 
distributions generally have luminosities 2-4 orders of magnitude 
lower (e. g. Biermann et al. 1981). Second, this explanation does 
not account for the morphology St the X-ray structure (e. g. , that the 
two arms are both pointed towards the nucleus of M87)» nor for its 
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spatial association with the radio halo (e. g. , that the position 
angles agree with the principal radio axes). There is no doubt that 
the radio halo is associated with M87, as it is f- n .r too bright and 
large to reside at the distance of the background cluster, 

V CONCLUSIONS 

The principal observational results of this study are the HRI 
residual map (Figure 1b), the VLA 1,4 GHz map (Figure 2b), and their 
superposition (Figure 4). The asymmetrical X-ray emission and the 
radio halo .are clearly associated with each other, but are not 
coincident. The ratio of X-ray to radio surface brightness ranges by 
a factor of 15, the structures differ in shape, and prominent radio 
features show no corresponding X-ray emission. Despite the quality 
of the observational material, the physical interpretation of the 
structures is not obvious. The radio halo clearly consists of 
energetic particles ejected from the nucleus. Its morphology is 
somewhat unusual, but it probably can be classified as a wide angle 
tail (WAT) source (Appendix B) . The principal issue concerns the 
cause of the associated X-ray emission. If inverse Compton is the 
operative mechanism then the average magnetic field must be at least 
10 times weaker than the equipartition or minimum energy value, and 
must vary considerably across the source. The alternative 
explanation is that the ejection of the radio emission plasma has 
caused local enhancement of thermal bremsstrahlung from the dense 
ambient medium. This should occur is the radio jet traveled outward 
at supersonic velocities, shocking and compressing the ambient 
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medium. The X-ray excess would persist for ~1 x 10 years before 
it cooled or dissipated. 

We can envision two experiments that might distinguish between 
the inverse Compton and thermal models, The most definitive test 
would be an X-ray spectrum of the asymmetrical component. Spatial 
resolutions < 30" and spectral resolutions E/A E > 100 could determine 
whether the excess emission is a power law or is comprised of 
emission lines. This observation is well within the capability of 
the planned Advanced X-ray Astrophysics Facility and perhaps other 
missions. A less definitive test, though one which could be 
accomplished today, would be narrow band optical imagery around the 
Ha line. This would determine whether cool gas is present in the 
vicinity of the east and southwest X-ray structures. Cool emission 
line filaments are known to exist in the densest gas within 40* of 
the nucleus (Ford and Butcher 1 9 7 9 > » but no detailed searches have 
been conducted in the halo region. While a negative result would not 
be informative, the detection of optical emission lines adjacent to 
the X-ray components would support a thermal model. 
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Appendix A 

Eccentricity of the X-ray Emitting Gas 

In their original analysis of the M87 HRX image, SGF reported 

that the X-ray surface brightness isophotes beyond the asymmetric 

component associated with the radio halo are approximately circular 

in contrast to the considerable eccentricity of the optical isophotes 

in this region of M87. However, after the subtraction of the 

symmetric X-ray component described in section Ha above, large-scale 

regions of low surface brightness emission appeared along a 

northwest-southeast axis, It is best seen in Figure 1c as 

enhancements at the level of the lowest contour. V band photographic 

surface photometry of the stellar emission in M87 shows an 

eccentricity of about 10% around 1' from the nucleus increasing to 

about 40% around 3* , and remaining at this high level beyond 30* 

(Carter and Dixon 1978). The orientation of the major axis is 

constant at P. A. 155°+ 5°, which is consistent with the X-ray 

enhancements. The X-ray excess appears in two regions, a ~8" 

diameter semicircle centered ~4' southeast of the nucleus, and a 5' x 

10' region centered ~7' northwest of the nucleus. (Recall that 10' 

corresponds to about M kpc. ) The mean surface brightness is a few 

percent of the symmetric component, and the total luminosity of the 

41 -1 

regions is ~5 x 10 ergs s 
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Ellipticity of X-ray isophotes In clusters or galaxies has not 
been widely discussed. Fabricant et al. (1984) have analyzed one 
case, Abell 2256, where the entire X-ray distribution (not just a few 
percent of the surface brightness) is elliptical with an axial ratio 
of 1,5 - 2.0 on scales of 3' to 30* <100 to 1000 kpc). The discovery 
of X-ray eccentricity in M87 is not surprising, as the X-ray gas must 
respond to the shape of the gravitational potential, which is 
presumably closely linked with the stellar distribution. A detailed 
analysis, which is outside the scope of this paper, of the X-ray and 
optical eccentricities might give insight into the distribution of 
non-luminous matter in that galaxy, The mass-to-luminosity ratio 10* 
from the center of M87 is >50, as determined from the gravitational 
confinement of the symmetric X-ray gas (Fabricant and Gorenstein 
1983), indicating that non-luminous mass dominates the gravitational 
potential at that radius. If the non-luminous mass has a less 
eccentric distribution than the stars, as may be suggested by the 
weak X-ray eccentricity, then a stellar origin of the non-luminous 
mass might be excluded, 

Appendix B 

Virgo A as a Radio Galaxy 

We examine here the VLA maps of the radio halo (Figures 2a-c ) 
not from the perspective of the associated X-ray emission, but as an 
improved view of an unusual radio galaxy. As described in section 
Ilia, the halo might be divided into four structures: (a) the knot 

2. 2' west of the nucleus; (b) the knot or lobe 3* -4' to the east; (c) 
the elongated component 2* -5' to the southwest; and (d) the large 
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scale envelope surrounding the entire souroe. He present here two 
rather different explanations of this str *cture, drawn from recent 
discussions of the origins of radio source morphologies. 

One interpretation involves ejection of inhomogeneous opposing 
jets by a precessing nucleus. Models of this type have been proposed 
for radio galaxies and quasars possessing S-shaped symmetry (see 
Ekers 1982), If one includes the inner lobes of the 'core' component 
of Vir A, the the principal features of the source can be Interpreted 
as three pairs of radio lobes ejected by a nucleus rotating counter 
clockwise. .The first ejection was along the NE/5W axis and produced 
halo components (c) and <d) which, being oldest, extend furthest from 
the nucleus and are most diffuse. The second ejection produced the 
components (a) and (b) along a E/W axis. This third ejection is 
still in progress, and is responsible for the radio/optical/X-ray 
synchrotron jet and f he inner lobes along an ESE/UNU axis. This 
explanation emphasizes the symmetries of the halo and would attribute 
asymmetries (e. g. , components c and d are not equally luminous; 
components a and b are not equidistant from the nucleus) to 
interactions between the ejected material and ambient medium. 

The second interpretation involves the continuous ejection of 
jets which bend or twist in various directions. Vir A is one of a 
number of radio galaxies (frequently cD giant ellipticals) at the 
center of rich clusters, and hence lie at rest in the center of a 
dense intracluster gas. The radio galaxies of these galaxies often 
show distorted structures that do not readily fit into standard 
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morphological classification schemes (classical doubles, head-tails, 
etc.). Examples of these central galaxies with radio extents similar 
to that of the Vir A halo <~60 kpc) include: 3C338 in Abell 2199, 
which has steep spectrum ridges and lobes (Burns et al. 1983); Per A 
= 3C8H in the Perseus cluster, which has an amorphous low surface 
brightness outer lobe surrounding a bright core and jet (Gisler and 
Miley 1 979); and 3 C 3 1 8 . 1 in MKW 3S, which has only one lobe with a 
very steep spectrum (Fomalont et al. in preparation). Though each is 
unique in detail, they have been generally classified as wide angle 
tail sources (WATs) characterized by ejecta that travel slowly 
outward from the central galaxy, bending and twisting in 
unpredictable directions (Burns 1983). Their morphologies might be 
determined by buoyancy forces, motions or inhomogeneities in the 
external medium, or electric currents traveling along their jets. If 
Vir A were in a distant cluster of galaxies and only a low resolution 
map like Figure 2a were available, it probably would be classified as 
a WAT with two distorted jets. The twisting jet model is supported 
by the apparent continuity of the surface brightness and the magnetic 
field orientation in Figures 2b and 3 that suggest a direct 
connection between the west knot (component a) and the southwest jet 
(component c). Such a connection would not be expected in the 
precessing nucleus model. 

Elucidation of the source physics in Vir A would benefit from 
high quality spectral index maps and a detailed multi-frequency 
polarization study to unravel the complex mixture of intrinsic 
polarized emission and foreground Faraday effects that are likely to 
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be present. While such data are not yet available, the polarization 
map shown in Figure 2c is very intriguing in this regard. The 
southwest component is strongly polarized, with polarization 
fractions up to 28%, while the eastern component is totally 
unpolarized, with an average polarization <0.8% <3 sigma limit). 

This polarization distribution is consistent with that found by 
Seielstad and Weiler <1969) in a one-dimensional scan and by 
Andernach et al. (1 979). but differs from that of Kotanyi (1 980). 
Andernach et al. report a polarization fraction of 21 + 5% for the 
south-west component at H. 8 GHz, with the peak polarization intensity 
appearing ~3' to the southwest of the peak total intensity. Kotanyi 
finds a peak polarization of 6 % for the east component at 610 MHz 
and significant polarized flux from the vicinity of the western knot 
{our component a), but little polarized from the southwest region. 
It is not clear whether these differences are instrumental in origin, 
result from different beam widths, or are intrinsic to the source. 

The reason for the startling difference in polarization fraction 
between the east and southwest, sides of the halo in our 1 . M GHz map 
is not self-evident. Perhaps the eastern component have detailed 
polarization substructure that is blurred by the 39" beam. 
Alternatively, it may be oriented away from us so that its radiation 
must pass through the dense gaseous core, subjecting its polarized 
flux to many small scale rotations that randomize the polarization 
angles. The southwestern structure would then be directed towards us 
and suffer less Faraday rotation. The apparent homogeneity of its 
field structure in Figures 2a and 2b suggests that there is not much 
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differential Faraday rotation across the component. Comparison of 
the position angle P. A. ~ 100°- 160° of the polarized flux in out 1.4 
GHz map (which is consistent with the 147°+ 14° angle found by 
Seielstad and Ueiler [19691) with the higher frequency polarization 
orientations reported by Andernach et al. (1979) indicates that 
considerable Faraday rotation is present. We believe it is 
premature to estimate the effect quantitatively, given the complex 
structure and incompatible beam sizes. 


Table 1 


Radio and X-ray Surface Brightness at 
Selected Locations in the Halo 


Region 

CJ$ 

Polar. 

<%> 

North radio halo 

0. 27 

<3 

East radio peak 
and X-ray spot 

1. 00 

<1 

East X-ray peak 

0. 65 

<1 

Southeast X-ray 
bridge 

0. 36 

1 

West radio .knot 

1. 07 

1 

Southwest X-ray 
peak 

0. 68 

6 

Southwest radio 
peak 

1 . 00 

5 

Southwest radio 
polarization peak 

0. 55 

28 


HRI 

cts/beam 

S x 

(yjy) 

log (S x 

/si. 4 : 

) a 

rx 

<24 

<0. 005 

<-7 • 

72 

V 

! 

O 

94 

81 

0. 017 

-7. 

76 

-0. 

94 

1 68 

0. 036 

-7. 

26 

-0. 

88 

75 

0. 01 6 

-7. 

35 

-0. 

89 

<24 

<0. 005 

< -8. 

31 

>-1 . 

01 

1 21 

0. 026 

~7. 

42 

-0. 

90 

44 

0. 009 

-8. 

03 

-0. 

97 

2A 

0. 005 

-8. 

02 

-0. 

97 
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Figure Captions 


Figure 1. Einstein HRI map of Vir A = M87: <a) the original image 

smoothed with a 15" FWHM Gaussian (contour levels 1, 2, 3, 4 , 6, 8, 

1 0, 15. 20, 50, and 90% of the peak 220 HRI counts/beam); <b> the 
residual emission after subtraction of the symmetrical component, 
smoothed to If*" FWHM (contour levels -0.75. 0.75, 1.5, 2.5, 3.5, 5, 

1 5, 50, 90% of the peak 218 HRI counts/beam); and (c) the same 
residual map smoothed to 39" FWHM for comparison with the radio maps 
(contour lev.els -1, 1, 2, 3, 4, 6, 10, 50, and 90% of the peak 2.61 x 
1 0 3 HRI counts/beam). 

Figure 2. VLA 1.4 GHz map of vir A = M87: (a) tapered map with 72" 

FWHM synthesized beam (contour levels -0.2, 0.2, 0.5, 1.0, 1.5, 2.0, 
2.5, 5, 1 0, 20, 30, 50, 70, and 90% of the peak 108 Jy/beam; 
polarization vector length 1* corresponds to 0.19 Jy/beam)i and (b) 

untapered map with 39" FWHM beam (contour levels -0.1, 0.1, 0.3, 0.3, 

0.6, 0.9, 1.2, 1.5, 5, 20, 50, and 90% of the peak 72.4 Jy/beamj 

polarization vector length 1' corresponds to 0.15 Jy/beam). 

Figure 3. VLA 5 GHz map of Vir A = M87, tapered with 40" FWHM 
synthesized beam (contour levels -0.2, -0.1, 0.1, 0.2, 0.3, 1, 5, 20, 

50, and 90% of the peak 42.5 Jy/beam). 

Figure 4. Superposition of selected contours of the 39" resolution 
radio (solid lines) and X-ray (dashed lines) maps, with the eight 
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locations listed in Table 'i marked. 
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